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(57)Abstract: 

PROBLEM TO BE SOLVED: To provide an air-fuel ratio 
control device capable of maintaining good controlling 
performance by properly setting the upper and lower 
limitations of the air-fuel ratio control amount for each 
cylinder in which the factoKs) likely to influence the 
accuracy of the by-cylinder air-fuel ratio presumed by 
the use of an observer is taken into consideration. 
SOLUTION: When the rotational variation amount MEMF 
of an engine is below the rotational variation threshold 
MFJUDOBS, control is made so that the allowable range 
decided by the upper and lower limitations KOBSVLMH 
and KOBSVLML of the by-cylinder corrective factor 
KOBSV#N widens gradually (S242, S249, S252). The 
upper limitation KOBSVLMH is controlled to be below 
the upper limit guard value KOBSLMHD while the lower 
limitation KOBSVLML is controlled to be over the lower 
limit guard value KOBSLMLD (S250, S251, S253, S254). 
wherein the upper limit guard value KOBSLMHD and the 
lower limit guard value KOBSLMLD are set in 
compliance with the atmospheric pressure PA. 
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* NOTICES * 

JPO and NCIPI are not responsible for any 
damages caused by the use of this translation. 

1 .This document has been translated by computer. So the translation may not reflect the original 
precisely. 

2.**** shows the word which can not be translated. 

3. In the drawings, any words are not translated. < 



CLAIMS 



(57) [Claim(s)] 

[Claim 1] OBUSABA which observes the internal state based on the model which describes the 
behavior of said engine's exhaust air system to be the air-fuel ratio detection means formed in an 
internal combustion engine's exhaust air system is set up. An air-fuel ratio presumption means 
classified by gas column to presume the air- fuel ratio of each gas column by considering the output 
of said air- fuel ratio detection means as an input, In the air-fiiel ratio control system of the internal 
combustion engine having an Air Fuel Ratio Control means classified by gas column to compute the 
amount classified by gas column of Air Fuel Ratio Control for carrying out feedback control of the 
air- fuel ratio of the gaseous mixture supplied to said each gas column so that the air- fuel ratio of 
each presumed this gas column may be completed as desired value The air-fiiel ratio control system 
of the internal combustion engine characterized by having a limit processing means to process so 
that it may enter in the tolerance where a bound value is set up according to either [ at least ] said 
engine's amount of rotation fluctuation, or atmospheric pressure, and said amount classified by gas 
column of Air Fuel Ratio Control is decided by said bound value. 



[Translation done.] 
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DETAILED DESCRIPTION 



[Detailed Description of the Invention] 
[0001] 

[Field of the Invention] Especially this invention relates to the air-fiiel ratio control system which 
carries out feedback control of the air- fuel ratio of the gaseous mixture supplied to an engine 
according to a gas column by the feedback control adapting the observer based on modern control 
theory about an internal combustion engine's air-fuel ratio control system. 
[0002] 

[Description of the Prior Art] It is prepared in an internal combustion engine's exhaust air system set 
section, and the air- fuel ratio presumption equipment classified by gas column which presumes the 
air- fuel ratio according to gas column using an observer is conventionally known based on the output 
of the air-fuel ratio sensor which generates the output proportional to an air- fuel ratio (for example, 
JP,6-173755,A). With this equipment, when the presumed air-fuel ratio becomes a value with a 
predetermined bound value out of range, it has prevented that a presumed air- fuel ratio emits by 
returning a presumed air-fuel ratio to initial value (value equivalent to A/F14.7). 
[0003] 

[Problem(s) to be Solved by the Invention] However, although emission of a presumed air- fuel ratio 
could be prevented by returning to initial value with the above-mentioned conventional equipment 
when a presumed air- fuel ratio became out of range [ a predetermined bound value ], there was a 
problem that converging again took time amount. Moreover, although the rate of flow of exhaust air 
affects presumption of the air-fuel ratio according to gas column, with the above-mentioned 
conventional equipment, the condition or air density (atmospheric pressure) of rotation fluctuation of 
the engine which affects the exhaust air rate of flow are not taken into consideration. 
[0004] Moreover, when feedback control of the air-fiiel ratio of the gaseous mixture which computes 
the amount of Air Fuel Ratio Control according to the detected air- fuel ratio, and is supplied to an 
engine using the amount of Air Fuel Ratio Control is carried out and the amount of Air Fuel Ratio 
Control exceeds a predetermined bound value, generally setting the amount of Air Fuel Ratio 
Control as the predetermined bound value is performed. However, the predetermined bound value in 
this case was a fixed value, and when applied to the feedback control of the air- fuel ratio according 
to gas column according to an air-fiiel ratio according to presumed gas column which described this 
technique above as it is, depending on the condition and air density (atmospheric pressure) of 
rotation fluctuation of an engine, said predetermined bound value turned into an unsuitable value, 
and it had the case where controllability ability was worsened. 

[0005] This invention is made in view of the point mentioned above, the bound value of the amount 
of Air Fuel Ratio Control according to gas column is more appropriately set up in consideration of 
the factor which affects the precision of the air- fuel ratio according to gas column presumed using an 
observer, and it aims at offering the air- fuel ratio control system which can maintain good 
controllability ability. 
[0006] 

[Means for Solving the Problem] An air- fuel ratio detection means by which this invention was 
prepared in an internal combustion engine's exhaust air system in order to attain the above- 
mentioned purpose, An air-fuel ratio presumption means classified by gas column to set up 
OBUSABA which observes the internal state based on the model which describes the behavior of 



http://www4.ipdl.ncipi.go.jp/cgi-bin/tran_web_cgi_ejje 



11/6/2006 



JP,3046948,B [DETAILED DESCRIPTION] 



Page 2 of 13 



said engine's exhaust air system, and to presume the air- fuel ratio of each gas column by considering 
the output of said air-fuel ratio detection means as an input, In the air- fuel ratio control system of the 
internal combustion engine having an Air Fuel Ratio Control means classified by gas column to 
compute the amount classified by gas column of Air Fuel Ratio Control for carrying out feedback 
control of the air- fuel ratio of the gaseous mixture supplied to said each gas column so that the air- 
fuel ratio of each presumed this gas column may be completed as desired value A bound value is set 
up according to either [ at least ] said engine's amount of rotation fluctuation, or atmospheric 
pressure, and it is characterized by having a limit processing means to process so that it may enter in 
the tolerance where said amount classified by gas column of Air Fuel Ratio Control is decided by 
said bound value. 

[0007] According to this configuration, based on the output of an air- fuel ratio detection means, the 
air-fuel ratio of each gas column is presumed. The amount classified by gas column of Air Fuel 
Ratio Control for carrying out feedback control of the air- fuel ratio of the gaseous mixture supplied 
to each gas column so that the air- fuel ratio of each presumed this gas column may be completed as 
desired value is computed. The amount of Air Fuel Ratio Control It is processed so that it may enter 
in the tolerance decided by the bound value set up according to either [ at least ] an engine's amount 
of rotation fluctuation, or atmospheric pressure. Consequently, the suitable bound value according to 
the presumed precision of the air- fuel ratio according to gas column can be set up, and good 
controllability ability can be maintained. 
[0008] 

[Embodiment of the Invention] The gestalt of operation of this invention is explained with reference 
to a drawing below. 

[0009] Drawing 1 is drawing showing the configuration of the internal combustion engine 
(henceforth an "engine") concerning one gestalt of operation of this invention, and its control unit. 
One is the engine of a 4-cylinder among this drawing. 

[0010] The inlet pipe 2 of an engine 1 is open for free passage to the combustion chamber of each 
gas column of an engine 1 through a tee (inlet manifold) 11. The throttle valve 3 is allotted in the 
middle of the inlet pipe 2. The sensor 4 is connected with the throttle valve 3 whenever [ throttle 
valve-opening ] (thetaTH), the electrical signal according to thetaTH is outputted whenever [ throttle 
valve-opening ], and the electronic control unit (henceforth "ECU") 5 is supplied. The supplementary 
air path 6 which bypasses a throttle valve 3 is established in the inlet pipe 2, and while being this 
path 6, the amount control valve 7 of supplementary airs is allotted. The amount control valve 7 of 
supplementary airs is connected to ECUS, and the amount of valve opening is controlled by ECU5. 
[001 1] The upstream of the throttle valve 3 of an inlet pipe 2 is equipped with the intake-air 
temperature (TA) sensor 8, and the detecting signal is supplied to ECU5. The chamber 9 is formed 
between the throttle valve 3 of an inlet pipe 2, and the inlet manifold 11, and the absolute-pressure 
(PBA) sensor 10 of inhalation of air is attached in the chamber 9. The detecting signal of the PBA 
sensor 10 is supplied to ECU5. 

[0012] The body of an engine 1 is equipped with the engine water temperature (TW) sensor 13, and 
the detecting signal is supplied to ECUS. The location sensor 14 is connected to ECU5 whenever 
[ crank angle / which detects angle of rotation of the crankshaft (not shown) of an engine 1 ], and the 
signal according to angle of rotation of a crankshaft is supplied to ECUS. The gas column distinction 
sensor by which the location sensor 14 outputs [ whenever / crank angle ] a signal pulse (henceforth 
a "CYL signal pulse") whenever [ predetermined crank angle / of the specific gas column of an 
engine 1 ] in a location, The top dead center at the time of charging-stroke initiation of each gas 
column It is related with (TDC). Whenever [ front crank angle ] in a location [ whenever / 
predetermined crank angle ] ((every) a 4-cylinder engine - 180 crank angles) from the CRK sensor 
which generates one pulse (henceforth a "CRK signal pulse") a fixed crank angle period (for 
example, periods of 30 degrees) shorter than the TDC sensor and TDC signal pulse which output a 
TDC signal pulse - changing - a CYL signal pulse - A TDC signal pulse and a CRK signal pulse 
are supplied to ECUS. These signal pulses are used for detection of various timing control, such as 
fuel injection timing and ignition timing, and an engine speed NE. 

[0013] The fuel injection valve 12 is formed in the upstream for every gas column for a while, each 
injection valve is electrically connected to ECU5 while connecting with the fuel pump of the inlet 
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valve of an inlet manifold 1 1 which is not illustrated, and fuel injection timing and fuel injection 
duration (valve-opening time amount) are controlled by the signal from ECU5. The ignition plug 
(not shown) of an engine 1 is also electrically connected to ECU5, and ignition timing thetalG is 
controlled by ECUS. 

[0014] The exhaust pipe 16 is connected to the combustion chamber of an engine 1 through the tee 
(exhaust manifold) 15. the part to which tees 15 gather to an exhaust pipe 16 — the broader-based 
air- fuel ratio sensor (henceforth a "LAF sensor") 1 7 is immediately formed in the downstream. 
Furthermore the directly under three way component catalyst 1 9 and the under floor three way 
component catalyst 20 are allotted to the downstream of the LAF sensor 17, and it is equipped with 
the oxygen density sensor (henceforth "02 sensor 1 ') 18 among these three way component catalysts 
19 and 20. Three way component catalysts 19 and 20 purify HC, CO, NOx, etc. in exhaust gas. 
[0015] It connects with ECU5 through the low pass filter 22, and the LAF sensor 17 outputs the 
electrical signal which carried out proportionally [ abbreviation ] to the oxygen density in exhaust 
gas (air- fuel ratio), and supplies the electrical signal to ECUS. 02 sensor 18 has the property that the 
output changes rapidly before and after theoretical air fuel ratio, and the output serves as a high level 
from theoretical air fuel ratio by the rich side, and it serves as a low by the Lean side. 02 sensor 18 is 
connected to ECU5 through the low pass filter 23, and the detecting signal is supplied to ECUS. 
[0016] The exhaust air reflux device 30 is established in the middle of the exhaust air reflux way 31 
which connects the chamber 9 and exhaust pipe 16 of an inlet pipe 2, and the exhaust air reflux way 
31, detects whenever [ valve-opening / of the exhaust air reflux valve (EGR valve) 32 which controls 
the amount of exhaust air reflux, and the EGR valve 32 ], and consists of the lift sensor 33 which 
supplies the detecting signal to ECUS. The EGR valve 32 is a solenoid valve which has a solenoid, 
and it connects with ECU5, and the solenoid is constituted so that whenever [ valve-opening ] can 
make it change with the control signals from ECUS. 

[0017] An engine 1 has at least the valve timing change-over device 60 of the high-speed valve 
timing which was suitable for the engine high-speed rotation field in the valve timing of an inlet 
valve, and the low-speed valve timing suitable for a low-speed rotation field switchable to two steps 
among an inlet valve and an exhaust valve. He stops one side of the two inlet valves at the time of 
low-speed valve timing selection, and is trying for a change-over of this valve timing to secure 
further the combustion stabilized when an air-fuel ratio was Lean-ized from theoretical air fuel ratio 
the time also including a change-over of the amount of valve lifts. 

[0018] The solenoid valve and oil pressure sensor (not shown) which the valve timing change-over 
device 60 performs a change-over of valve timing through oil pressure, and perform this oil pressure 
change-over are connected ECUS. The detecting signal of an oil pressure sensor is supplied to 
ECUS, and ECU5 controls a solenoid valve and performs change-over control of valve timing. 
[0019] Moreover, the atmospheric pressure (PA) sensor 21 which detects atmospheric pressure is 
connected to ECUS, and the detecting signal is supplied to ECU5. 

[0020] ECUS operated orthopedically the input signal wave from the various sensors mentioned 
above, corrected a voltage level to predetermined level, and is equipped with the input circuit which 
has the function of changing an analog signal value to a digital signal value, a central-process circuit 
(a CPU), the store circuit which consists of ROM and RAM which memorizes various maps 
[ program / operations / various / and ] mentioned later, the result of an operation, etc. performed by 
this CPU, and the output circuit which output a driving signal to the various solenoid valves and the 
ignition plug of fuel injection valve 12 grade. 

[0021] According to an engine operation condition, ECUS calculates the fuel injection duration 
TOUT of a fuel injection valve 12 with the following formula 1, and outputs the signal which drives 
a fuel injection valve 12 based on this result of an operation while it distinguishes various engine 
operation conditions according to the output of the LAF sensor 17 and 02 sensor 18, such as a 
feedback control operating range and an opening control operating range, based on various above- 
mentioned engine operation parameter signals. 
[0022] 

[Equation 1] TOUT=TIMFxKTOT ALxKCMDMxKL AFxKOB S V#N drawing 2 is a functional 
block diagram for explaining the calculation technique of the fuel injection duration TOUT with the 
above-mentioned formula 1, and explains the outline of the calculation technique of the fuel 
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injection duration TOUT in the gestalt of this operation with reference to this. In addition, although 
the fuel amount of supply to an engine is computed as fuel injection duration with the gestalt of this 
operation, since this corresponds to the fuel quantity injected, it is also calling TOUT fuel oil 
consumption or fuel quantity. 

[0023] In drawing 2 , block Bl computes the basic fuel quantity TIMF corresponding to an 
inhalation air content. Although this basic fuel quantity TIMF is fundamentally set up according to 
an engine speed NE and the absolute pressure PBA of inhalation of air, it is desirable to perform 
amendment which modeled the inhalation-of-air system from a throttle valve 3 to the combustion 
chamber of an engine 1 , and took the delay of inhalation air into consideration based on that 
inhalation-of-air system model. In that case, thetaTH and atmospheric pressure PA are further used 
whenever [ throttle valve-opening ] as a detection parameter. 

[0024] Block B-2 - B4 are multiplication blocks, carry out the multiplication of the input parameter 
of a block, and output it. By these blocks, the operation of the above-mentioned formula 1 is 
performed and fuel oil consumption TOUT is obtained. 

[0025] By carrying out the multiplication of all the feedforward system correction factors, such as 
the purge correction factor KPUG set up according to purge fuel quantity at the time of the purge 
activation by the EGR correction factor KEGR set up according to the amount of exhaust air reflux 
during the engine water temperature correction factor KTW set up according to the engine water 
temperature TW, and exhaust air reflux activation, and the evaporation fuel processing unit, block 
B9 computes a correction factor KTOTAL, and inputs it into block B-2. 

[0026] Block B21 determines the target air-fiiel ratio multiplier KCMD according to engine-speed 
NE, the absolute pressure PBA of inhalation of air, etc., and inputs it into block 22. Since the target 
air-fiiel ratio multiplier KCMD is, are proportional to the inverse number, i.e., fiiel-air-ratio F/A, of 
air- fuel ratio A/F, and a value 1 .0 is taken at the time of theoretical air fuel ratio, it is also called 
target equivalent ratio. Block B22 corrects the target air- fuel ratio multiplier KCMD based on 02 
sensor output VM02 inputted through a low pass filter 23, and inputs it into blocks B18 and B23. 
Block B23 performs fuel cooling correction according to a KCMD value, computes the policy 
objective air-fuel ratio multiplier KCMDM, and inputs it into block B3. 

[0027] Block B10 samples the LAF sensor output value inputted through a low pass filter 22 for 
every generating of a CRK signal pulse, carries out the sequential storage of the sampled value at 
ring buffer memory, chooses the sampled value sampled to the optimal timing according to the 
engine operation condition (LAF sensor output selection processing), and it inputs it into block B18 
through the low pass filter block B16 while inputting into block Bl 1. This LAF sensor output 
selection processing takes into consideration that the air-fuel ratio which changes depending on the 
timing of a sampling is correctly undetectable, and that time amount until the exhaust gas discharged 
from a combustion chamber reaches the LAF sensor 17, and the reaction time of the LAF sensor 
itself change with engine operation conditions. 

[0028] Block Bl 1 has a function as the so-called observer, presumes the air- fuel ratio for every gas 
column based on the air-fiiel ratio of the set section (mixed gas of the exhaust gas discharged from 
each gas column) detected by the LAF sensor 17, and inputs it into the blocks B12-B15 
corresponding to four gas columns. In drawing 2 , block B12 corresponds to gas column #1, block 
B13 corresponds to gas column #2, block B 14 corresponds to gas column #3, and block B15 
corresponds to gas column #4. The air-fuel ratio (air- fuel ratio which the observer block B12 
presumed) of each gas column computes correction factor KOBSV#N classified by gas column 
(N=l-4) by PID control so that it may be in agreement with a set section air- fuel ratio, and blocks 
B12-B15 input it into block B5-B8, respectively. 

[0029] Block Bl 8 computes the PID correction factor KLAF by PID control according to the 
deflection of a detection air- fuel ratio and a target air- fuel ratio, and inputs it into block B4. 
[0030] As mentioned above, with the gestalt of this operation, while applying the PID correction 
factor KLAF which the output of the LAF sensor 17 responded and was computed by the usual PID 
control to the above-mentioned formula 1, the fuel oil consumption TOUT for every gas column (N) 
is computed by applying further correction factor KOBSV#N classified by gas column set up 
according to the air-fiiel ratio of each gas column presumed based on the LAF sensor output to the 
above-mentioned formula 1 . Dispersion in the air- fuel ratio for every gas column can be canceled by 
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correction factor KOBSV#N classified by gas column, the rate of purification of a catalyst can be 
raised, and a good exhaust gas property can be acquired in various engine operation conditions. 
[0031] With the gestalt of this operation, since data processing by CPU of ECUS realizes, the 
function of each block of drawing 2 mentioned above explains the contents of processing concretely 
with reference to the flow chart of this processing. In addition, in the following explanation, a suffix 
(k) corresponds to the sampling time in a discrete time system, and (k) etc. corresponds to a value 
etc. a value and last time this time, respectively (k-1). However, especially (k) that shows a value 
this time is omitted as long as there is no need. 

[0032] Drawing 3 is the flow chart of the processing which computes the PID correction factor 
KLAF and the correction factor KOBSV according to gas column according to the output of the LAF 
sensor 17. This processing is performed for every generating of a TDC signal pulse. 
[0033] At step SI, it distinguishes whether it is [ ****** ] starting mode, i.e., under [ cranking ], and 
shifts to processing (step S10) of starting mode at the time of starting mode, if it is not starting mode 
- calculation (step S2) of the target air- fuel ratio multiplier (target equivalent ratio) KCMD and the 
policy objective air-fuel ratio multiplier KCMDM, and LAF sensor output selection processing — 
carrying out (step S3) — detection equivalent ratio KACT is calculated (step S4). The detection 
equivalent ratio KACT changes the output of the LAF sensor 17 into equivalent ratio. 
[0034] Subsequently, activity distinction of whether activation of the LAF sensor 17 was completed 
is performed (step S5). As compared with a predetermined value (for example, 0.4V), this 
distinguishes the difference of the output voltage and the main electrical potential difference of the 
LAF sensor 17 as activation was completed when this difference was smaller than a predetermined 
value. 

[0035] Next, an engine operation condition distinguishes whether it is in the operating range 
(henceforth a "LAF feedback field") which performs feedback control based on the output of the 
LAF sensor 17 (step S6). When an engine operation condition is in a LAF feedback field, it is set as 
"1", and the feedback control flag FLAFFB is set as "0", when other. And air- fuel ratio correction 
factor KOBSV#N classified by gas column and the PID correction factor KLAF are calculated (steps 
S7 and S8), and this processing is ended. When it is feedback control flag FLAFFB=1, at step S7, 
correction factor KOBSV#N classified by gas column is computed by the processing (processing 
after step S203) of drawing 9 mentioned later, and by step S8, the PID correction factor KLAF is 
computed so that the detection equivalent ratio KACT may be in agreement with the target 
equivalent ratio KCMD with well-known PID control. When it is feedback control flag FLAFFB=0, 
in step S7, correction factor KOBSV#N classified by gas column is set as "1.0", and is set as the 
predetermined value determined according to an engine operation condition at step S8. 
[0036] Drawing 4 is the flow chart of the LAF feedback field distinction processing in step S6 of 
drawing 3 . 

[0037] At step SI 21, it distinguishes first whether the LAF sensor 17 is in an inactive condition. The 
flag FFC which shows by "1" that it is [ fuel ] under cut when it is in an active state distinguishes 
whether it is "1" (step SI 22), and when it is FFC=0 The flag FWOT which shows by "1" that it is 
[ throttle-valve ] under full open distinguishes whether it is "1" (step SI 23), and when it is not 
FWOT=l Distinguish (step SI 24), and when it is VBAT>=VBLOW, whether battery voltage VBAT 
detected by the sensor which is not illustrated is lower than the predetermined lower limit VBLOW 
It distinguishes whether there is any gap (LAF census TOIKI gap) of the LAF sensor output 
corresponding to theoretical air fuel ratio. And when one answer of steps S 121 -S 125 is affirmation 
(YES), the feedback control flag FLAFFB which shows by "1" that feedback control based on a LAF 
sensor output can be performed is set as "0" (step SI 32). 

[0038] On the other hand, when all of the answer of steps S121-S125 are negation (NO), it judges 
with activation of the feedback control based on a LAF sensor output being possible, and the 
feedback control flag FLAFFB is set as " 1 " (step S 1 3 1 ). 

[0039] Next, calculation processing of correction factor KOBSV#N classified by gas column in step 
S9 of drawing 3 is explained. 

[0040] The presumed technique of the air- fuel ratio according to gas column by the observer is 
explained to the beginning, and the calculation technique of correction factor KOBSV#N classified 
by gas column according to the air-fiiel ratio according to gas column presumed next is explained. 
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[0041] It considered that the air-fiiel ratio of the exhaust air system set section was the weighted 
average in consideration of the time contribution of the air- fuel ratio of each gas column, and the 
value at the time of time of day k was expressed like a formula 2. In addition, fuel quantity (F) is 
written as a control input, and fuel-air-ratio F/A is used with the formula 2. 
[0042] 
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That is, although the fuel-air ratio of the set section carried out the multiplication of weight C (40% 
and a that front are 30%, — , etc. for the gas column which burned immediately before) to the 
combustion hysteresis of the past for every gas column, it was expressed in total. If a block diagram 
expresses this model, it will become like drawing 5 and that equation of state will become like a 
formula 3. 
0043] 

^Equation 3] 



x (k-2)1 fOlOlfx (k-3)1 |"0l 
x (k-l) = 0 0 1 I I x (k- 2) + 0 u 
x (k) J Lo 0 OJ L* Ck- 1) J L 1 J 



(k) 



Moreover, if the fuel-air ratio of the set section is set with y (k), the method of an output can express 
an equation like a formula 4. 
[0044] 
[Equation 4] 



y (k) = [c 



fx (k-3) 1 
iCacJ x (k-2) +C4 u(k) 
L* (k-l) J 



ci : 0. 05, c a : O. 15, c 8 : 0 . 30, 
c* : 0 . 0 5 t-?Z>. 

In a formula 4, since it cannot observe, even if u (k) designs an observer from this equation of state, 
x (k) cannot observe it. Then, if [ assuming that the air- fuel ratio in front of 4TDC (namely, the same 
gas column) is in the steady operation condition of not changing rapidly ] x(k+l) = x (k-3), a formula 
4 will become like a formula 5. 
[0045] 
[Equation 5] 
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x (k-2) 
* (k-1) 

x (k) 
x (k + 1 ) 



0 10 0 
0 0 10 
0 0 0 1 
10 0 0 



x (k-3) 
x (k-2) 
x (k - 1 ) 
x (k) 



y (k) = [cic jc 8 c 4 ] 



(k-3) 
(k-2) 
(k-1) 
x (k) 



Thus, it is checked experimentally that the set-up model models the exhaust air system of a 4- 
cylinder engine well. Therefore, the problem which presumes the air-fuel ratio according to gas 
column comes back from set section A/F to the problem of the usual Kalman filter that the equation 
of state shown with a formula 6 and the method of an output observe x (k) by the formula. If the load 
matrices Q and R are set like a formula 7 and Riccati's equation is solved, the gain matrix K will 
become like a formula 8. 
[0046] 
[Equation 6] 

X (k + 1 ) =AX (k) +B u (k) 
y (k) =CX (k) +Du .(k) 



A = 



0 


1 


0 0 


0 


0 


1 0 


0 


0 


0 1 


.1 


0 


0 0 



C = [cic 2 C3C 4 ] B=D= [0] 



X (k) = 




[0047] 
[Equation 7] 
f 1 oo ol 
Q J o ioo 

W 0 0 10 1 
|_0 0 0 1 J 



R= [1] 



[0048] 
[Equation 8] 

- 0.3 0 9 3 
1 . 19 16 
0. 3 0 9 3 
O . 0 8 0 3 



K = 



In the model of this operation gestalt, since there is no input u (k) in the configuration of a general 
observer, as shown in drawing 6 , it becomes the configuration which considers y (k) as an input, and 
if this is expressed with a formula, it will become like a formula 9. 
[0049] 
[Equation 9] 
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X (k + l) = CA-KC] X (k) +Ky (k) 

-AX (k) + K (y (k) -CX <k) ) 

x (k-3) " 
£ (k-2) 

4 ck-i) 
£ (k) 

Therefore, the estimate X hat (k) of this fuel-air ratio according to gas column is computable from 
the estimate X hat (k) of set section fuel-air-ratio y (k) and the past fuel-air ratio according to gas 
column. 

[0050] When computing the fuel-air-ratio X hat classified by gas column (k+l) using the above- 
mentioned formula 9, although the detection equivalent ratio KACT (k) is applied, CX hat (k) and 
the (weighting aggregate value of the fuel-air ratio according to [ of four ] gas column) do not 
include delay as set section fuel-air-ratio y (k) to this detection equivalent ratio KACT (k) including 
the response delay of the LAF sensor 17. Therefore, in having used the formula 9, it is the effect of 
the response delay of the LAF sensor 17, and the fuel-air ratio according to gas column cannot be 
presumed correctly. When especially the engine speed NE is high, since the recurrence interval of a 
TDC signal pulse becomes short, the effect of response delay becomes large. 

[0051] So, with this operation gestalt, the estimate X hat (k+l) of the fuel-air ratio according to gas 
column was computed by computing the estimate y hat (k) of a set section fuel-air ratio with a 
formula 10, and applying this to a formula 11. 
[0052] 

[Equation 10] 

(k) =DLy (k-l> + (1-DL) CX (k) 

[0053] 

[Equation 11] 

£ <k+ 1) =AX (k) +K (y (k) -y* <k) ) 

In the above-mentioned formula 10, DL is a parameter equivalent to the time constant of the 
response delay of the LAF sensor 17, and is computed using DL table shown in drawing 7 with this 
operation gestalt. DL table is set up so that DL value may turn into a value between 1.0 from 0 
according to an engine speed NE and the absolute pressure PBA of inhalation of air. In this drawing, 
PBA(s) 1-3 are 660mmHg 5 460mmHg, and 260mmHg, and they compute the detected engine speed 
NE and the time constant DL according to the absolute pressure PBA of inhalation of air by 
performing a interpolation operation suitably, for example, respectively. In addition, it is checked 
experimentally that the value of a time constant DL has the optimal value equivalent to time amount 
later about 20% than the value equivalent to an actual response time delay. Moreover, since change 
of as opposed to change of the absolute pressure PBA of inhalation of air in a time constant DL is 
small, you may make it choose the average value of the value of the time constant DL corresponding 
to PBA1, PBA2, and PBA3 only according to an engine speed NE, for example so that clearly from 
drawing 7 . 

[0054] In addition, in formulas 10 and 1 1, all the values of a component (x hats (k-3), x hats (k-2), x 
hats (k-1), x hats (k)) make the initial vector of X hat (k) the vector of the detection equivalent ratio 
KACT, and also let initial value of y hat (k-1) be the detection equivalent ratio KACT in a formula 
10. 

[0055] Thus, by using the formula 1 1 which transposed CX hat (k) in a formula 9 to the estimate y 
hat (k) of a set section fuel-air ratio including the response delay of a LAF sensor, the response delay 
of a LAF sensor can be compensated appropriately and the exact air- fuel ratio according to gas 
column can be presumed. In addition, presumed equivalent ratio KACT#l(k) -KACT#4(k) of each 
gas column in the explanation which is the following is equivalent to x hats (k), respectively. 
[0056] Next, the technique of computing correction factor KOBSV#N classified by gas column 
based on the presumed air- fuel ratio according to gas column is explained with reference to drawing 
8. 
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[0057] First, as shown in a formula 12, the desired value KCMDOBSV as equivalent ratio 
corresponding to target A/F (k) is computed by doing the division of the detection equivalent ratio 
KACT corresponding to set section A/F with the last operation value of the average of correction 
factor KOBSV#N classified by gas column of all gas columns. # The gas column exception of a 1 
cylinder, calculate correction factor KOBSV#l by PID control so that deflection DKACT#l(k) 
(=KACT#l(k)-KCMDOBSV (k)) with the desired value KCMDOBSV (presumed equivalent ratio 
KACT#l(k of k) and # 1 cylinder) may be set to 0. 
[0058] 

[Equation 12] 

KCMDOBSV (k) =— j KACT (k) 

£ KOB S V#N /4 

N- 1 

Proportional KPOB#l, integral term KIOB#l, and differential term KDOB#l are calculated with a 
formula 13, and, more specifically, correction factor KOBSV#l is further computed a gas column 
exception with a formula 14. 
[0059] 

[Equation 13] KPOB#l(k) =KPOBSVxDKACT#l (k) 
KIOB#l(k) =KIOBSVxDKACT#l(k)+KIOBS#l (k-1) 
KDOB#l(k) =KDOBSVx (DKACT#l(k)-DKACT#l (k-1)) 
[0060] 

[Equation 14] KOBSV#l(k)=KPOB#l(k)+KIOB#l(k)+KDOB#l(k) 

# Perform an operation with the same said of 2 - # 4-cylinder, and compute KOBS V#2-#4. 

[0061] Thereby, since the air- fuel ratio of each gas column is converged on a set section air-fuel ratio 

and a set section air- fuel ratio is converged on a target air-fuel ratio with the PID correction factor 

KLAF, the air- fuel ratio of all gas columns can be completed as a target air- fuel ratio as a result. 

[0062] Furthermore, correction factor study value KOBSVR#N classified by gas column which is 

the study value of this correction factor KOBSV#N classified by gas column is memorized to RAM 

which computed for every operating range and was backed up by the following formula with the dc- 

battery. 

[0063] 

[Equation 15] KOBSVR#j=CRxKOBSV#N+(l-CR) xKOBSVR#Nj - KOBSVR#Nj of the 
operating-range parameter (j=l-3) determined that j will mention later here according to the absolute 
pressure PBA of inhalation of air, the weighting factor with which CR is set as the value between 0 
and 1 , and the right-hand side is the last study value. 

[0064] Drawing 9 is the flow chart of the correction factor KOBSV#N calculation processing 
classified by gas column in step S7 of drawing 3 (Maine). 

[0065] First, at step S201, the feedback control flag FLAFFB distinguishes whether it is "1", it is 
FLAFFB=0, and when there is nothing to a LAF feedback field, correction factor KOBSV#N 
classified by gas column is set as "1.0" (step S202), initialization processing of step S226 is 
performed, and this processing is ended, while setting all (k-1) of x hats (k-3) of said formulas 10 
and 1 1, x hats (k-2), x hats (k-1), x hats (k), and y hat of a formula 10 as the detection equivalent 
ratio KACT at step S226 - DKACT#N of a formula 13 - "0" moreover, initialization processing 
which sets KIOB#N of a formula 13 as correction factor KOBSV#N classified by gas column is 
performed. 

[0066] When it is FLAFFB=1 at step S201, a operating range is distinguished according to the 
absolute pressure PBA of inhalation of air (refer to steps S203-S207 and drawing 14 ). Namely, 
when it is the low loading field where the absolute pressure PBA of inhalation of air is lower than 
the 1st place constant pressure PBOBRF 1 The operating-range parameter j is set as "1" (steps S203 
and S205). The absolute pressure PBA of inhalation of air by the one or more 1st place constant 
pressures PBOBRF And when it is an inside load field lower than the 2nd place constant pressure 
PBOBRF 2 higher than the 1st place constant pressure PBOBRF 1 The operating-range parameter j 
is set as "2" (steps S203, S204, and S206), and when the absolute pressure PBA of inhalation of air is 
a two or more 2nd place constant pressures [ PBORF ] heavy load field, the operating-range 
parameter j is set as "3" (steps S203, S204, and S207). 
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[0067] At continuing step S208, KOBSVLMH/L decision processing ( drawing 10 ) is performed 
and the upper limit KOBSVLMH and lower limit KOBSVLML of correction factor KOBSV#N 
classified by gas column are determined. And an engine speed NE is the predetermined bound values 
NOBSVH and NOBSVL (for example, 4000 rpm, respectively). ** is distinguished [ whether 
whether it being within the limits of lOOOrpm and the absolute pressure PBA of inhalation of air are 
within the limits of the predetermined bound values PBOBSVH and PBOBSVL (for example, 
respectively 760mmHg(s), 200mmHg), and ] (refer to step S209 and drawing 14 ). NE<=NOBSVL, 
Or when it is NE>=NOBSVH, PBA<=PBOBSVL, or PBA>-PBOBSVH Correction factor 
KOBSV#N classified by gas column is set as study value KOBSVR#Nj corresponding to a operating 
range at that time (step S221), and limit processing of correction factor KOBSV#N classified by gas 
column is performed (steps S222-S225). That is, when it considers as KOBSV#N=KOBSVLMH 
when correction factor KOBSV#N classified by gas column is larger than a upper limit 
KOBSVLMH (steps S222 and S225), it considers as KOBSV#N=KOBSVLML when smaller than a 
lower limit KOBSVLMLH (steps S223 and S224), and it is within the limits of a bound value, it 
progresses to said step S226 as it is. 

[0068] When the answer of step S209 is affirmation (YES), it distinguishes whether the feedback 
control flag FLAFFB was "1" at the time of the last activation of this processing (step S210), and 
when it is last FLAFFB=0, it progresses to said step S221, and initializing is performed. When it is 
FLAFFB=1 also last time, KOBSV#N calculation processing (factice) of drawing 1 1 and 
KOBSVR#N calculation processing of drawing 12 are performed (steps S21 1 and S212), and this 
processing is ended. 

[0069] Drawing 10 is the flow chart of the KOBSVLMH/L decision processing in step S208 of 
drawing 9 . 

[0070] At step S241, the rotation fluctuation threshold MFJUDOBS is computed with the following 

formula 16. 

[0071] 

[Equation 16] 

MFJUDOBS=MFJUDxKMFJDOBS - here, MFJUD is the basic value of the rotation fluctuation 
threshold computed by searching the MFJUD map set up according to an engine speed NE and the 
absolute pressure PBA of inhalation of air, and the MFJUD map is set [, so that an engine speed NE 
increases, and ] up so that the absolute pressure PBA of inhalation of air increases, and the basic 
value MFJUD may decrease. Moreover, KMFJDOBS is the correction factor of a rotation fluctuation 
threshold, and this correction factor KMFJDOBS is computed by searching the KMFJDOBS table 
set up according to the absolute pressure PBA of inhalation of air. The KMFJDOBS table is set [ as 
shown in drawing 15 (a), ] up so that the absolute pressure PBA of inhalation of air increases, and a 
correction factor KMFJDOBS may decrease. 

[0072] At continuing step S242, it distinguishes whether the detected amount MEMF of engine 
rotation fluctuation is larger than the threshold MFJUDOBS computed at step S241. When this time 
value of the time amount (namely, amount proportional to 1/NE) taken for an engine crankshaft to 
carry out fixed include-angle (for example, 30 degrees) rotation here is set to ME (k), the amount 
MEMF of rotation fluctuation is a parameter defined by the following formula. 
[0073] MEMF=ME(k)-ME(k- 1 ) 

When the amount MEMF of rotation fluctuation is larger than a threshold MFJUDOBS It 
distinguishes whether correction factor KOBSV#N classified by gas column of the gas column 
(henceforth a "this time gas column") set as the object of an operation this time is smaller than the 
Lean side predetermined value KOBSVMFL (for example, 0.96) (step S243). When it is 
KOBSV#N<KOBSVMFL and the value of correction factor KOBSV#N classified by gas column 
inclines toward the Lean side All (j=l-3) of a upper limit KOBSVLMH, a lower limit KOBSVLML, 
and study value KOBSVR#Nj are set as "1.0" (step S244). Predetermined time TOBSVMF (for 
example, 3 minutes) is made to set and start the down count timer tmOBSVMF (step S245), and this 
processing is ended. 

[0074] By steps S242-S244, when correction factor KOBSV#N classified by gas column of a gas 
column inclines toward the Lean side greatly and this time, engine rotation fluctuation Since all of 
the bound values KBOSVLMH and KOBSVLML and study value KOBSVR#Nj are set to "1.0" 



http://www4.ipdl.ncipi.go.jp/cgi-bin/tran_web_cgi_ejje 11/6/2006 



JP,3046948,B [DETAILED DESCRIPTION] 



Page 11 of 13 



Correction factor KOBSV#N classified by gas column is also set as "1.0" (non-correction value) by 
processing of steps S278-S283 of drawing 1 1 mentioned later, and the increment in the rotation 
fluctuation by the air- fuel ratio of a specific gas column being Lean-ized too much can be prevented. 
[0075] The answer of steps S242 or S243 is negation (NO), i.e., MEMF<=MF JUDOB S , and even if 
it is amount MEMF>MF JUDOB S of rotation fluctuation when the amount of rotation fluctuation is 
small or, the value of the timer tmOBSVMF started at step S245 when correction factor KOBSV#N 
classified by gas column was larger than the Lean side predetermined value KOBSVMFL 
distinguishes whether it is "0." And tmOBSVMF> While being 0, this processing is ended 
immediately and it is set to tmOBSVMF=0, and when a negative (NO) condition continues [ the 
answer of steps S242 or S243 ] beyond the predetermined time TOBSVMF, the value of the down 
count timer tmDKOBSV started at the following step S248 distinguishes whether it is "0" (step 
S247). tmDKOBSV> While being 0, when this processing is ended immediately and it is set to 
tmDKOBSV=0, predetermined time TDKOBSV (for example, 10 seconds) is made to set and start 
Timer tmDKOBSV (step S248). 

[0076] Subsequently, a upper limit KOBSVLMH is updated by adding the predetermined value 
DKOBSVMF (for example, 0.001) (step S249), and it distinguishes whether it is smaller than the 
upper limit guard value KOBSLMHD to which the this updated upper limit KOBSVLMH is set 
according to atmospheric pressure PA (step S250). And when it is KOB S VLMH<KOB S LMHD , 
when it is KOB S VLMH>=KOB SLMHD again, immediately, a upper limit KOBSVLMH is set as 
the upper limit guard value KOBSLMHD (step S251), and it progresses to step S252. 
[0077] At step S252, a lower limit KOBSVLML is updated by subtracting the predetermined value 
DKOBSVMF, and it distinguishes whether it is larger than the minimum guard value KOBSLMLD 
to which the this updated lower limit KOBSVLML is set according to atmospheric pressure PA (step 

5253) . And when it is KOBSVLML>KOBSLMLD, when it is KOBSVLML<=KOBSLMLD again, 
immediately, a lower limit KOBSVLML is set as the minimum guard value KOBSLMLD (step 

5254) , and this processing is ended. 

[0078] Here, the upper limit guard value KOBSLMHD and the minimum guard value KOBSLMLD 
are determined in the guard value calculation processing (processing by which processing which 
synchronized with the TDC signal pulse is not performed and which is performed in the so-called 
background) shown in drawing 13 by searching the guard value table shown in drawing 15 (b) 
according to an atmospheric pressure PA. The guard value table is set up so that the upper limit 
guard value KOBSLMHD may decrease, so that atmospheric pressure PA falls, and the minimum 
guard value KOBSLMLD may increase, namely, so that the tolerance of the bound values 
KOBSVLMH and KOBSVLML may become narrow, in drawing 15 (b), the predetermined values 
KOBSLMHD 1 and KOBSLMHD2 are set, for example as 1.02 and 1.05, respectively - having - 
the predetermined values KOBSLMLD 1 and KOBSLMLD2 - respectively - 0. - it is set as 98 and 
0.95. 

[0079] Step S By 248 or less processing, the increment of a upper limit KOBSVLMH and the 
decrement of a lower limit KOBSVLML are performed for every predetermined time TDKOBSV, 
and the tolerance of correction factor KOBSV#N classified by gas column is expanded gradually. 
And if the amount MEMF of rotation fluctuation exceeds a threshold MFJUDOBS and correction 
factor KOBSV#N classified by gas column of a gas column becomes smaller than the Lean side 
predetermined value KOBSVMFL by processing of steps S242-S244 this time, both the bound 
values KOBSVLMH and KOBSVLML will be returned to "1.0." 

[0080] Drawing 1 1 is the flow chart of the KOBSV#N calculation processing (factice) in step S21 1 
of drawing 9 . 

[0081] First, at step S271, the air- fuel ratio according to gas column is presumed by the technique 
mentioned above (calculation of presumed equivalent ratio KACT#N), and, subsequently 
proportional KPOB#N, integral term KIOB#N, and differential term KDOB#N of the feedback 
control of air- fuel ratio according to gas column are computed with said formula 13 (step S272). 
Limit processing of integral term KIOB#N is performed at continuing steps S273-S276. That is, 
when it is set as the upper limit KOBSVLMH when integral term KIOB#N is larger than a upper 
limit KOBSVLMH, it progresses to step (steps S273 and S274) S283, it is set as the lower limit 
KOBSVLML when integral term KIOB#N is smaller than a lower limit KOBSVLML, and it 
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progresses to step (steps S275 and S276) S281 and integral term KIOB#N is within the limits of the 
bound values KOBSVLMH and KOBSVLML, it progresses to step S277. 
[0082] In step S277, correction factor KOBSV#N classified by gas column of a gas column is 
computed with said formula 14 this time, and limit processing of correction factor KOBSV#N 
classified by gas column is performed at continuing steps S278-S283. Namely, when correction 
factor KOBSV#N classified by 1 gas column is larger than a upper limit KOBSVLMH integral term 
KIOB#N - last time - value KIOB#N (k-1) returning (step S282) - Correction factor KOBSV#N 
classified by gas column is set as a upper limit KOBSLMH (step S283). When correction factor 
KOBSV#N classified by 2 gas columns is smaller than a lower limit KOBSVLML integral term 
KIOB#N - last time -- value KIOB#N (k-1) - returning (step S280) - When correction factor 
KOBSV#N classified by gas column is set as a upper limit KOBSVLML (step S281) and correction 
factor KOBSV#N classified by 3 gas columns is within the limits of the bound values KOBSVLMH 
and KOBSVLML, this processing is ended as it is. 

[0083] While correction factor KOBSV#N classified by gas column is computed according to 
presumed equivalent ratio KACT#N classified by gas column according to processing of drawing 
11 , limit processing of correction factor KOBSV#N classified by gas column and integral term 
KIOB#N which were computed is carried out so that it may enter within the limits of the bound 
values KOBSVLMH and KOBSVLML. 

[0084] Drawing 12 is the flow chart of the correction factor study value KOBSVR#Nj calculation 
processing classified by gas column in step S212 of drawing 9 . 

[0085] In this drawing, at step S401, it distinguishes whether predetermined time progress was 
carried out from this feedback control initiation point in time classified by gas column (at the shift 
time [ The execution condition of the feedback control according to gas column ] to the condition 
[ condition / abortive ] of execution condition formation), this processing is immediately ended 
before progress, and after progress progresses to step S402. At step S402, equivalent ratio deflection 
DKACTOB#N which is the deflection of the detection equivalent ratio KACT (k) and equivalent 
ratio KACT#N classified by presumed gas column (k) is computed by the following type. 
[0086] DKACTOB#N=KACT(k)-KACT#N(k) 

And this processing is ended, without absolute value |DKACTOB#N| of equivalent ratio deflection 
distinguishing whether it is below the predetermined deflection DKOBRFLM (step S403), and 
computing study value KOBSVR#Nj, when it is |DKACTOB#N|>DKOBRFLM. 
[0087] Moreover, when it is |DKACTOB#N|<=DKOBRFLM, study value KOBSVR#j is computed 
for every operating range as follows, and this processing is ended. 

[0088] That is, when the absolute pressure PBA of inhalation of air is lower than the 1st place 
constant pressure PBOBRF 1, study value KOBSVR#Nl is computed with said formula 15, and it 
memorizes to Backup RAM (steps S404 and S406). Moreover, the absolute pressure PBA of 
inhalation of air is the one or more 1st place constant pressures PBOBRF, and when lower than the 
2nd place constant pressure PBOBRF 2, study value KOBSVR#N2 is computed with said formula 
15, and it memorizes to Backup RAM (steps S404, S405, and S407). Moreover, when the absolute 
pressure PBA of inhalation of air is the two or more 2nd place constant pressures PBOBSRF, study 
value KOBSVR#N3 is computed and it memorizes to Backup RAM (steps S404, S405, and S408). 
[0089] As explained in full detail above, with this operation gestalt, the engine amount MEMF of 
rotation fluctuation in the condition below the predetermined threshold JFJUDOBS While the 
tolerance decided by the bound values KOBSVLMH and KOBSVLML of correction factor 
KOBSV#N classified by gas column is expanded gradually If the amount MEMF of rotation 
fluctuation exceeds a threshold MFJUDOBS and correction factor KOBSV#N classified by gas 
column of a gas column becomes smaller than the Lean side predetermined value KOBSVMFL this 
time Since both the bound values KOBSVLMH and KOBSVLML were returned to "1.0", suitable 
tolerance (bound value) according to the engine amount MEMF of rotation fluctuation can be set up. 
Moreover, since the guard values KOBSLMHD and KOBSLMLD of the bound values 
KOBSVLMH and KOBSVLML were set up according to atmospheric pressure PA, a suitable bound 
value can be set up according to the altitude the car with which the engine 1 was carried runs. 
Consequently, change of the exhaust air rate of flow which affects the precision of the presumed air- 
fuel ratio according to gas column is taken into consideration, and the stable control can be 
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maintained by narrowing tolerance at the time of the increment in the amount of rotation fluctuation 
with the inclination for presumed precision to fall, or high ground. 

[0090] With this operation gestalt, step S271 of drawing 11 is equivalent to the air- fuel ratio 
presumption means classified by gas column, correction factor KOBSV#N classified by gas column 
is equivalent to the amount classified by gas column of Air Fuel Ratio Control, steps S272 and S277 
of drawing 1 1 are equivalent to the Air Fuel Ratio Control means classified by gas column, and 
processing of drawing 10 and drawing 13 is equivalent to a limit processing means at steps S278, 
S279, and S281 of this drawing, and S283 list. 

[0091] In addition, do not restrict this invention to the operation gestalt mentioned above, and the 
guard values KOBSLMHD and KOBSLMLD of the bound values KOBSVLMH and KOBSVLML 
are good also as a fixed value. Or the ** which does not perform expansion processing (steps S247- 
S254) of the tolerance decided by the bound values KOBSVLMH and KOBSVLML according to the 
amount MEMF of rotation fluctuation, You may make it set up the bound value KOBSVLMH and 
the KOBSVLML itself like the inclination shown in drawing 15 (b) according to atmospheric 
pressure PA. 
[0092] 

[Effect of the Invention] As explained in full detail above, according to this invention, based on the 
output of an air-fuel ratio detection means, the air- fuel ratio of each gas column is presumed. The 
amount classified by gas column of Air Fuel Ratio Control for carrying out feedback control of the 
air- fuel ratio of the gaseous mixture supplied to each gas column so that the air- fuel ratio of each 
presumed this gas column may be completed as desired value is computed. The amount of Air Fuel 
Ratio Control Since it is processed so that it may enter in the tolerance decided by the bound value 
set up according to either [ at least ] an engine's amount of rotation fluctuation, or atmospheric 
pressure, the suitable bound value according to the presumed precision of the air- fuel ratio according 
to gas column can be set up, and good controllability ability can be maintained. 



[Translation done.] 
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* NOTICES * 

JPO and NCIPI are not responsible for any 
damages caused by the use of this translation. 

l.This document has been translated by computer. So the translation may not reflect the original 
precisely. 

2 **** s hows the word which can not be translated. 
3. In the drawings, any words are not translated. 



DESCRIPTION OF DRAWINGS 



[Brief Description of the Drawings] 

[Drawing 1] It is drawing showing the configuration of the internal combustion engine concerning 1 
operation gestalt of this invention, and its air- fuel ratio control system. 

[Drawing 2] It is a functional block diagram for explaining the Air Fuel Ratio Control technique in 
this operation gestalt. 

[Drawing 3] It is the flow chart of the processing which computes an air-fuel ratio correction factor 
based on a LAF sensor output. 

[Drawing 4] It is the flow chart of LAF feedback field distinction processing. 

[Drawing 5] It is the block diagram of the model in which the behavior of an internal combustion 

engine's exhaust air system is shown. 

[Drawing 6] It is the block diagram showing the configuration of the observer in this operation 
gestalt. 

[Drawing 7] It is drawing showing the table for setting up the response delay time constant (DL) of a 
LAF sensor. 

[Drawing 8] It is a block diagram for explaining the feedback control of air- fuel ratio according to 
gas column. 

[Drawing 9] It is the flow chart of the main routine which computes the correction factor 
(KOBSV#N) according to gas column. 

[Drawing 10] It is the flow chart of the processing which determines the bound value of the 
correction factor according to gas column. 

[Drawing 1 1] It is the flow chart of the subroutine which computes the correction factor according to 
gas column. 

[Drawing 12] It is the flow chart of the processing which computes the study value of the correction 
factor according to gas column. 

[Drawing 13] It is the flow chart of the processing which computes the guard value used by 
processing of drawing 10 . 

[Drawing 14] It is drawing showing the operating range which performs feedback control of air-fiiel 
ratio according to gas column. 

[Drawing 15] It is drawing showing the table used by processing of drawing 10 or drawing 13 . 
[Description of Notations] 

1 Internal Combustion Engine (Body) 

2 Inlet Pipe 

5 Electronic Control Unit (ECU) (Air-fuel Ratio Presumption Means Classified by Gas Column, Air 
Fuel Ratio Control Means Classified by Gas Column, Limit Processing Means) 
12 Fuel Injection Valve 

16 Exhaust Pipe 

17 Broader-based Air-fiiel Ratio Sensor (Air- fuel Ratio Detection Means) 
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«ncD£^(i 'Oj *LT. MISfS'J3S*SS 

JtMiEfS&KOBS V?XS'JP I DfiiEf^StKLAF 
WiS^^ffoT (XT77S7. SS) . *3ail£i&7 
\ t rty9fflffl77>fFL AF FB= 1 T'« 
£££, X^-y^S 7-Cti^jd>-r-5ia9<7DS!lS Uf7 
ys 2 0 3^(DHgm) IZ&V. %ffiwmiE&%i.K O B 
SV#N*»J|Ca*n. Xr 7 7"S8T(t J^»<73PID 
$iJfP(CJ;f3, ^mSSJtKACT^B^SfiJtKCMD 
lC-Sfr-5J;5ICP I D«E^S:KLAF^*ttJ^tl 

7-r-hVt-yi7$iJ^77^FLAFFB=0T?»^> 
£^T«, Xf 7 7"S7T(l ^fSgiJWlE^ScKOBS V 
#N«» ri. 0j lCi£5££*U 7f 7 7'S 8TttI>y 
>3HEttfi8 fcfc C ^ n ^> fc $ tl ^ . 

[0 0 3 6] 04li, iSCXT'^S 6 fc*5^SLA 
F X ^ - hVt >y * ftWNgijgl&O 7 D-^^ - h 
•5. 

[0 0 3 7] 56fXf 7 7'S12im LAF-tr>+»- 

1 7&*m&vimz$>2>fr&fr&mmv, m&vtmz& 



( 5 ) ttlMft 3 0 4 6 9 4 8f 

9 10 

< 7 >f - h* A* v 9 Vm SJIfTWfil i #J5S L T, 7-f — H 
A*y^affll77yFLAFFB* r 1 j fCigJtf-S (X 
f7^S 13 1). 

[0 0 3 9] &»cH3©X'r.y:/S9fc*t**»«giJ*| 
iE^KOBS V#N«?Hffij8lgfC-3I,iTiftWr-5. 
[0 0 4 0] *!0l=3f7*-!f- i ''?»;«t*a«giJffi*81tO« 

C&ftff PltAjEfllftK O B S V # N©J*tH^ft*IJl9rr 

[0 0 4 1 ] *S5**^«©^*BJt*ftatt©^«ifc© 

k<B<t:S0>II*ft3!C2©.fc3l;:3iLfc. 
(F) **f£«i:l,fcfci&. KSC2Ttt«£J±F/A« 

[0 0 4 2] 

[ft 2] 

[F/A] (k) =C, [F/A# 1 ] +C 2 [F/A# 3] 

+ Ca [ F /A # 4 ] + C« [F/A# 2] 

[F/A] (k + 1) =C, [F/A» 3] +Ci [F/A#4] 
+C S [F/A#2] +C 4 [F/A#l] 

[F/A] (k+2) =C, [F/A#4] + C 8 [F/A»2] 
+ C S [F/A # 1] +C 4 [F/A#3] 

Km^C (0]*fcmifflK«5&L£:JSfS5fi4 0 %. [0 0 4 3] 

#3 0%. **J)tLfct>©0&ffl-T3SLfc. C [S3] 

co ; &7 ; ;U$'7'n-y H 5 ©J: -5 13^*0. 30 

x (k - 2 ) "J roiOlfx [k-3)"| ro"I 
x ( k - 1 ) = 0 0 1 x [ k - 2 ) + 0 u ( k ) 
x (k) J LoOOJLx (k-l)J L 1 J 

*£«<0*£tfc£y (k) ttS777jS:£ [0 0 4 4] 

fx (k-3) "I 
y (k) = C c t c )C|] I x (k-2) + c « u (k) 
|_x (k- 1) J 

c i : 0 . 05, c 2 : 0 . 15, c 8 : 0 . 30, 
c* : 0 . 0 5 k-?Z>. 

ftsuiz^x, u (k) \tmm^mx'$>^rctt>. e 
©«fi8*s^e.*7 r if-n<&is^LT'bx (k) \tm 

litC&-5£<K5tLTx (k+1) =x (k-3) fTS 
t, KSUttftiCSO.fc'SfcttS. 
[0-0 4 5] 

[ft 5 ] 50 



5fFFC*« r l j A»S^«fl9Jb Uf77"S 1 2 
2) . FFC=0TSI»t§ll XO-vHl/^HtT 
Wit* r l j TSt7 7^FWOTAi r 1 j 
*«S'JL (Xf7yS12 3) . FWOT=lT&lr>i 

VBAT#0f«TEIIVBLOW e fcO(g<->#S*^*IJSiJ 
L Ut77"S12 4) . VBAT^VBLOWlr^-5 

LAFt y^miKD-rti 
(LAFt>fxK*fn) atfcsa^^fcwsij-r 10 

-5. ^IT. Xf77"S 1 2 1~S 1 2 5(DU"rn^CD 
(YES) rotate. LAFt>tUIAl:i^ 

it7^-HA'7i'iW77?'FLAFFB* r 0 j \Z 
^-TS UfyT'S 13 2). 

[0 0 3 8] -7?- Xf77 , S12 1~S 1 2 5 
t^T§^ (NO) <7)<hi*«, L AF-fe>-9"tii73{-So* 



x (k-2> 
x (k - 1 ) 
x (k) 



_x (k + 1 ) J [_1 0 0 OJ _ 



11 

0 10 0* 
0 0 10 
0 0 0 1 



(k-3) 
(k-2) 
(k-1) 
x (k) 



y (k) = [cic 2 c 3 c 4 ] 



(k-3) 
(k-2) 
(k-1) 
x (k) 



( 6 ) 



:3 0 4 6 9 4 8^ 



12 



10 



(k) &m9t?zmiis<D*)U'?>7 4)u?<Dma\zMm 

ft*. 
[0 0 4 6] 
[ft6] 



X (k + 1 ) =AX (k) + B u (k) 
y (k) =CX (k) + Du (k) 



A = 



0 10 0 
0 0 10 
0 0 0 1 
10 0 0 



C = [ctc 2 c 3 c 4 ] B = D= [0] 



X (k) = 



x (k-3) 
x (k-2) 
x (k - 1 ) 
x (k) 



[0 0 4 7] 
[ft 7 ] 



[0 0 4 8] 
[ft 8] 



[1 00 0] 
0100 
0 0 10 
0 0 0 1 J 



R= [1] 



30 



f- 0. 3 0 9 3 
v _ I 1 . 19 16 
* 0. 3 0 9 3 

L 0. 0 8 0 3 



(k) ©**A*fr*«*tfto. cn*saT*t 

£ft5t9<£><fc?t;:ft£o 
[0 0 4 9] 
[ft 9] 

X Ck + 1) = [A-KC] X (k> + Ky (k) 
/N 

-AX (k) +K (y (k) -CX (k) ) 



40 



X (k) = 



x (k-3) 
$ (k-2) 
4 (k-1) 
£ (k) 



SJto*««XAyh (k) ^0<0»«*J*R£tfc 
©atjettx/wh (k) *#tbf 

[0 0 5 0] ±IEftit9^ffl^T^SiJ«S^ltXA^ h 
(k+1) *#ttrr*»^ m^SfB^^lty (k) £L 
T, tfctBSfiJtKACT (k) 3&*5iffl*n*^ £C0l& 
aSlltKACT (k) tt, LAFt>+M 7 

tts^T^scofcfcfu cxA^b (k) (4^<om 

CtttTSftv*. #^X>>?>IieftNE^ifi!^^# 
tt, TDC«^;i/X<0**IWIW*tfi<3a:*<0TJiS*iI 

n<z>ig»a**^<fts 0 

[0 0 5 1] fCT^lMTil ftjfcl OfcckO* 
£«JIS£it<0#£ttyAy h (k) S*ttiU ^MS 

sci 1 fcisfflf sc£kj;d, %mm^\t(Dm^my. 
/\vh (k+i) «rjrm-r*«fc"3fcbfc. 

[0 0 5 2] 
[ft 1 0] 

(k) =DLy (k-1) + (1-DL) CX (k) 

[0 0 5 3] 
[ftl 1] 



*^«*£Jty (k) &V»£(0«fff9Jj|B 50 



13 

X U+ 1) =AX (k) + K (y 



(k) 



-y (k) ) 



±fEftitl 0\Z&^T. DLtiL AFir>*y- 1 7 COJC^ 

y;Utt, DL«36^x>>?>|iIlEftNE&Of8ft»ert 
*»ffPBAl:«i;T0*5 1. O<Dra<0ffi<t&£,fc3 
llSSSnT^*. HHJw45^T, PBAl-3H*tl 
t'nMAH 6 6 0mmHg, 4 6 0mmHg, 2 60 
mmHgtiD, 3Sfil«ra{Rff*ff ^T, ft[lib/tx> 
v >B»6ft N E RtX»a«rt*6*tffi P B A KJg DfcWPJfe 
ftDLflDjnHSfTS. &*3, P*«*DL<0tttt, lUESCD 

PBA1, P B A 2 Rtf P B A 3 ^)P$St D L 

[0 0 5 4] tji& % ft5£l OS^l 1 fc*^T, XA7 

h (k) bMt. mzumigatx (x/wh 

(k-3) , x/Vyh (k-2) , x/\7h (k- 
1) , x/Vyh (k) ) ©tl#*T*lH3«ttKACT 
(B^HUtU S«10(:*^TyAyh (k-1) 

KCMDOBSV (k) =— — 



( 7 ) WflFfS 3 0 4 6 9 4 8^ 

14 

o»JHttt>«kaia«JtKACTi-r«. 

CO 0 5 5] COiiKl, ft^C9 KWt-5CXAy h 
<k) LAF*>iJ-0)j*:aF»tiS^ie*^««ffi 
lt©«JE«yAyh (k) K«#«*fc*SCl 
fcZltfcJ;?). LAF-fc>1*-0*aFiin*ia«^*i«L 
TiE«a««8iJffl«Jt©««*R5c:t38«T**. 

KT©R«k:*W4&*A«©Ji3eS*JtKACT# 
1 (k) -KACT#4 (k) -tft-^tlx/Vy h 
(k) fcWSf*. 
10 [0 0 5 6] *^li«Lfcai»9Jffij|(Rlt»C»'Cf^Tatt 
»J*|jE«»KOBSV#N*iHH'r**ttt, ^8£# 

[0 0 5 7] 5feT, 1 2 KS-T «t 5 *£rg&A/ 

F tC*HS-T *tftfflS«Jt K A C T &£Affi0>Aff»Jft|jE 
®KOBSV# N© IttlilttlFttTRJE ITS 

*A/Fl:»K«Silti:LT01«ltKCMDOB 
SV (k) *»tiJU # l««0)a«8UliiE«ftKOB 
SV#1H ^SMKCMDOBSV (k) t # 1 
ft«<B*£3*itKACT# 1 (k) t^IiDKAC 
20 T#l (k) ( = KACT#1 (k) — KCMDOBS 
V (k) ) j^OtWiil:, P I Dft9«K:J;0#«> 

[0 0 5 8] 

Ml 2] 
K AC T <k) 



£ KOB S V#N (k 

N- 1 



1) /4 



«»SK I O B # 1 Stf «»SK D O B # 1 $ 

4tckO»fB8iJ»iE«*KOBSV# 1*W 30 

[0 0 5 9] 

[ft 13] KPOB#l (k) =KPOBSVXDKA 
CT# 1 (k) 

KIOB#l (k) =K I OB S VXDKACT# 1 

(k) +KI OBS# 1 (k- 1) 
KDOB#l (k) =KDOBSVX (D K A C T # 1 

(k) — D K A C T # 1 (k - 1 ) ) 

[0 0 6 0] 

[ftl4]KOBSV#l (k)=KPOB#l (k) 40 
+ KIOB#l (k) +KDOB#l (k) 

# 2-u4^mz-D^x^>mm<Dmmn^, kobs 

V#2-#4£gLffit*£o 

[0 0 6 1] c*ik:J:0, ««tt0&*Jtl;M^ff&NI 
jtfcJR#U *-&«S!*8JtttP I DiiEfllRKLAFI: 

n<D^mit^BS^mit\ziRM^it^^^^^o 

[0 0 6 2] 2<b\Z, £<Z>afK»J«:E«ftKOBSV# 
N©^a«T*«««»J*liE«**SltKO B S V R # 



TAy^Ty^SftfcRAMCBttS. 
[0 0 6 3] 

115] KOBS VR# j =CRXKOBS V#X + 
(1-CR) XKOBSVR#Nj 
CdT% j ttttSB-rScti^Ka^lfiJtJEPBAHlE 

VTfcfetsnzmfcmw^? ( j = i - 3 ) , c 

R«O^bl0F B 1(OlH^$tll)i^m #iS<0K 
OBSVR#N j ttflftlllO^BfiT**. 
[0 0 6 4] H9lt HS^r-y^S 7 K is -5 SUB 

suffliE^ftKOB s vxNwmmm (*<» <D7u- 

[0 0 6 5] ftfXr7^S2 0 ITU 7^f-KA*7 
?1»77^FLAFFB^ r l j ^5^*JpJgiJ U F 
LAFFB = 0T»oTLAF7-<-hVty9««i:& 

»«9J*IjE«»KOBSV#N* ri. 0J 
izWifeL Uf7^S 2 0 2) , 7r>^S 2 2 60ft 

6TH ffiGftAl OXtXl ICO, xA^h (k - 

3) , x/Vyh (k-2) , x/Vyh (k-1) , xA 

y h (k) fttfRSCl 0<Dy/\y h (k-1) ^r^:T^ 

aa*jtKACT^K«-r*ttt>^, ssi3^dk 

ACT#N£ TOJ \Z. £fcfti£ 1 3 CDK I O B # N£ 
SNHttjEfllRKO B S V # N fCig^TS^JSIffc^LS* 



( 8 ) 

15 

[0 0 6 6] Xf7^/S 2 0 1TFLAFFB = 1 T & 
iRa«rtlft»jEEPBAfc*i;T31ft«*(0« 
SU^fT 1 ^ Ufy^S 2 0 3-S 2 0 7, @14# 
If) o -rftto-fe, »»«rtlft»fl£PBA)W*l»f3fejEP 
BOBRF 1 J:0«^«<MKr«l«T**i#»4, Slfcffi 

r l j KKJtl, Uf7/S 2 0 3, 
S 2 0 5) > ®ftf«J£PBA^i^j£PBO 
BRF1H±"C, *oil0f^EPBOBRFUOiS 
t«2BfSffPBOBRF 2 <k 0C^ + ft«f««T** 10 

r 2j ^isjeu (x^ 

7^S 2 0 3, S 2 0 4, S206) , Bft»«rt*J#flE 
P B ArfW, 2 BfSJE P B O R F 2 JK±<&iSfiffi««T* 

3HE««/X^^-^ j S r 3 J KKje-rs 
(Xf7^S 2 0 3, S 2 0 4, S 2 0 7 )o 
[0 0 6 7] a<Xf7^S 2 0 8T*H KOBSVL 

MH/L&^sa:® (0io) zmnis, §K.m»MJE%& 

KOBS V#N<D±K«SKOB S VLMHStfTElK 
OBSVLML^m o ^Lt, X>v>IeIte»N 
E#f3f«±TIB«lNOB S VH, NOBSVL (0ijx.tf 20 
•?-nfn4 0 0 0 r pm, 1 0 0 0 r pm) (D^mfalZ 
***>5*v &^©mi : PM*6*fffiPBA^0f^±TKM 
PBOBSVH, PBOBSVL (CTAtf-etl^n 7 6 
OmmHg, 2 0 0mmHg) (^ilrtC^^^?: 
SrWB'JL Ut7^S 2 0 9, H14#HS) , NE^N 
OBSVL, ifc«NE^NOBSVH, SfcliPBA 
^PBOBSVL, SitttPBA^PBOBSVHT* 
*<fc*tt, ««8"J*IjE«»KOBSV#N*» ^cd£# 

oaefi[«K:j*«-r**H«KOBsvR#Nj kr« 

L Uf7^S2 2 1) , gm'WjEMKOBSV# 30 
NCOU ^ >y h*QS£fT5 Uf7^S 2 2 2-S2 2 
5) . a«BiJ»jE«*KOBSV#N#±IB 
•(KOBSVLMHiO^Sl^Sft KOBSV#N 
= KOB S V L MH <h L (X^y 7S 2 2 2 , S22 
5) , TBUKOBSVLMLHiO/hSlitSH, K 
OBS V#N=KOBS VLMLchl (Xf7^S2 2 
3, S 2 2 4) , ±TRttO«lrtl:*5iSH t<D 
^SWIBX^^yS 2 2 6lCjttTo 

[0 0 6 8] Xf7^S 2 0 9W*^ (YES) <0 

*^OM(Hl*ff^^^^T^^-hV\*^^$fJ 40 
«77^FLAFFB* T»o £*JgiJ I, 

Uf7^S2 10) > MIH FLAFFB = 0T^o/: 
£*HU (HEXrv^S 2 2 1 fci^Ttt«MM85g*fr 
ffi@fcFLAFFB=lT»o&<h*tt, Ell 1© 
KOBS V#N^ffifi!l31 RZSmi 2<DKOBS 

VR#N»Wiaa*^fTbT Ut'^S 2 1 1, S2 
12), *Sa31^^T-r^)o 

[0 0 6 9] B10H 09(DXf7^S 2 0 8 K*5tt 
§KOB S VLMH/Lifc£fflS<B:7n-^-- h~C$> 

50 



4#WFSS 3 0 4 6 9 4 8^ 
16 

[0 0 7 0] Xf7^S2 4 1T1t TIB*j$16KJ; 
0, IsMg^l&HffiMF J UDOBS^gtilt^ 
[0 0 7 1] 
[»16] 

MFJUDOBS=MFJ U D X KM F J DOBS 

MFJUDU X>>*>0teftNERtf»lSff 
rt«S*tffiPBAfCJCi;TSftffiSnfcMF J UD7 7 7$ 

&tK MFJUD77-/H X>v>EKEftXE***JP 
-T *fc«a«rttt»JEPBA*«littl'r-S(a<!f, 
S*IMF JUD*tj^/>f ^ct^fCiS^^tlTC^o £ 

fc, kmfjdobsh imteSiftKfi^WiE^&T* 

0, £<Z>«IEfi*ScKMF J DOBSfct, KS*l^tft*fJE 
P BAKlS^TlS^tl/tKMF J DOB 

&m-rz>cti\z&owm£tiz>o kmfjdobst- 

0 15 (a) \Z^T£o\Z, »S<Bfrt*ft*NEP 
BA3W«lHI-r*«f*|jE«»KMF J DOBSW>f 

[0 0 7 2] g<Xf7^S 2 4 2 Til ftffil/:X> 
^>0(E3C»*MEMF*^ Xf7^S2 4 1T»HjL 
fcHfiMF J UDOBS£0*#V^S^£«S'J*r£o 
X>;x><£^>^Htft— (0iJA^3O 

tSi) O^lElffl^ME (k) IsMg^i&SM 
EMF«Ti3a^i$n5/l7^-^T**§o 
[0 0 7 3] M E M F =ME (k) -ME (k-1) 
@68»iMEMF^HfiMF J UDOBSi^Sli 

£^5) <^«fSBU*iE«»:KOB S V±r xa*'J — >ftl»r 
f IKOB SVMFL (OT;Ui'0. 9 6) «k9'J*Sl»a* 
5*^*WaiJL (Xf77S 2 4 3J . KOB5V = X< 
KOBS VMF LT*oTa«fgiJt*iE»aKOB S V =r 
NOfi^U->ffi|(*ffloTl^i:*!i, _h[gfgKOB S 
VLMFL TSI K O B S V L M L RL^KI KOBS 
VR#X J ( j =1-3) ££T r i. 0, KftifeL 
(XT77S 2 4 4) , ^>^7>^<-7tmOB 
S VMFSBfJtRrMTOBS VMF («A«3#) Kir 
yhLTX^-hili: (Xf7^S 2 4 5) , 

[0 0 7 4] Xf7yS 2 4 2-S244CJ:0, X> 

KOBSV#NWU->«|i:fi^T^5t*H ±T« 
IKBOSVLMH, KOBSVLMLR^fllKO 

BsvR#Nj*^t ri. oj tznzw. aje-r 

£0 1 lOXr^^S 2 7 8~S 2 8 3 (OfflSfclcfc DSl 
«g«J»iEffi»CKOBSV#Nt) Ti. 0j O&lftjEM) 

\zmmzn, &mf&n<D&mikibmmzv~>ik2ti2> 

[0 0 7 5] XT7^S 2 4 2 *fcttS 2 4 3 (Dgfj&tg 
£ (NO) , tftbSMEMFSMF J UDOB ST» 



( 9 
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oT@te£ih*a*/h;£^i£, *;fctt[lHE£ll>SMEM 
F>MF J UDOB ST*oTt>»l(f9J*iiE«*KOB 
S V # N ifi U - >fiij ffi m fit K O B S V M F L £ 0 * * ^ 
t#tt, Xfy^S 2 4 5TX^M/W7tmO 
BSVMFCDffl^ roj fr&fr&mq-fZo fLTtm 
OBSVMF>0T*§f B 1ll tt^Cttl^Tl, 
tmOBSVMF = 0 tftOXfy^S 2 4 2 SfcteS 
2 4 30*^Sfi (NO) ©ttaaWftWIBTOBSV 
MFHliJIlfcifft *W7f7yS 2 4 8T'X^ 
-h$n^>^^>^«>> h*-*f? tmDKOBS VOftt 10 
roj Jfc5*>**«j-r* (Xr'^S 2 4 7) . tm 

DKOBsv>oT$ > ^rait m^>\z^mm^mTL, 

tmDKOBSV=0tft§ < !:, ?^7tmDKOBS 
V^HHTDKOBSV m^iiom Ct7h 
LTX^-h£tl;£ Uf7yS 2 4 8) o 

[0 0 7 6] &V>-e, iMKOBSVLMH^, Bf^ 
MDKOB SVMF (OTAtfO. 0 0 1) *iraD£-r*c 
ii:«k0Mfb (XT7-/S 2 4 9) % »M«fLfc±RB 
IKOBSVLMH^, ^JBEPAKJ&^TlS^Sn* 
±(E^- FtKOBSLMHDctO/h^^S^^fflij 20 

(Xf7yS 2 5 0) o flTKOB SVLMH< 
KOBSLMHDT»5i*BE%t:, £fcKOBSV 
LMH^KOBSLMHDTWttlt ±ISfitKOB 
SVLMH^t^lK^- FiKOBSLMHDCtS 
IT Uf7^S2 5 1), Xf7^S 2 5 2 tCjttTo 

[0 0 7 7] Xf7^S 2 5 2T(1 TPSfitKOBSV 
LML£, ^lDKOBSVMF^»t§Cif:J; 
B»WU KH«fbfcTIH«KOBSVLML^ 
JEEP AKJKUTK^SftSTKtf-- KfitKOB S LML 
D±9;fc*^rtvg;fr**lJ8iJ-rs (Xr77 J S 2 5 3) = 30 
*LTKOBSVLML>KOBSLMLD"C*Si# 
tiB^fC SfcKOB S VLML^KOB SLMLDTf 
*5tSH T(SMKOBSVLML£^TPH#-K 
iKOBSLMLDCfSStt (XT77S 2 5 4) , 

[0 0 7 8] ddT, -tPfi^f — HIKOBSLMHDS 
^TB^-HttKOBSLMLDH m 1 3 \Z^Ttf- 

S) tC&^T, *&JEPAK:jSi;TH 1 5 (b) \Zxk~? 40 

-Hffl^-^«, *«BEPA**fiT"r-SHi:±liB^- 
FlKOBSLMHD»>L> Ti^-HIKOBS 
LMLD^*f}jQT£<£5K:, T&:b*>±TPgffiKOBS 
' VLMH, KOBSVLML©»$«l*«»<&-5ct5 
KtB^StiT^*. HI 5 (b) \Zte^Ts Sf^lKO 
BSLMHD1WKOBSLMHD2H -tn-etlW 
0 2, 1. 0 SfcltKJtSfrL BfatffiKOBSL 
MLD 1MKOBSLMLD2H ^rtl^tim^i 
0.9 8,0.9 5JCi8SSn-5. 50 



#fFS 3 0 4 6 9 4 8^ 
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[0 0 7 9] Xt7^S 2 4 86tT««ia^J:0, Bfjg 
^ffflTDKOBSV»:±SfiKOBSVLMH^ > 
2? U^>hWTMKOBSVLML©r^ U h 
***fTS*l, m^S'JttlE^^KOBS V#N(DfF#ftffl 
fit, Xf7^S 2 4 2-S2 
4 4(D5[iaiI^0, EieSliiMEMF^HlMFJU 
DOBS Sr«AA^^igaf8f©»«giJ»jE«ftKOB S 
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